Scaling of ground state fidelity in the thermodynamic limit: XY model and beyond 
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We study ground state fidelity defined as tlie overlap between two ground states of the same 
quantum system obtained for slightly different values of the parameters of its Hamiltonian. We 
focus on the thermodynamic regime of the XY model and the neighborhood of its critical points. We 
describe in detail cases when fidelity is dominated by the universal contribution reflecting quantum 
criticality of the phase transition. We show that proximity to the multicritical point leads to 
anomalous scaling of fidelity. We also discuss fidelity in a regime characterized by pronounced 
oscillations resulting from the change of either the system size or the parameters of the Hamiltonian. 
Moreover, we show when fidelity is dominated by non-universal contributions, study fidelity in the 
extended Ising model, and illustrate how our results provide additional insight into dynamics of 
quantum phase transitions. Special attention is put to studies of fidelity from the momentum space 
perspective. All our main results are obtained analytically. They are in excellent agreement with 
numerics. 



I. INTRODUCTION 

Quantum phase transitions (QPTs) happen when dra- 
matic changes in the ground state properties of a quan- 
tum system can be induced by a tiny variation of an ex- 
ternal parameter [1] . This system-specific parameter can 
be a magnetic field in spin systems 0, Q , intensity of a 
laser beam in cold atom simulators of Hubbard-like mod- 
els Q , dopant concentration in high-Tc superconductors 
[|, etc. 

Traditional condensed matter approaches to QPT fo- 
cus on identification of the order parameter and the pat- 
tern of symmetry breaking as well as on studies of two 
point correlation functions and the excitation gap jl| . At 
the critical point of the second order QPT the correlation 
length diverges while the gap in the excitation spectrum 
vanishes. This is typically described by power-law singu- 
larities. Correlation length ^ diverges as |A — Ad""^, while 
the excitation gap closes as |A — Acl^"^, where A is the ex- 
ternal field driving the transition, Ac marks the quantum 
critical point, and z and v are the critical exponents as- 
sociated with the universality class of the system. The 
exponents are universal in the sense that they do not 
dependent on the microscopic details of the system. 

A somewhat different ways of looking at QPTs have 
recently emerged from quantum information community. 
One of them is based on studies of quantum entanglement 
of spins/atoms/etc. undergoing a QPT M. Another is 
based on studies of ground state fidelity [3] (see [8] for 
a recent review). Finally, an extreme approach emerged 
where the Hamiltonian is designed to have a particular 
many-body state as its ground state .9]. We will study 
ground state fidelity below. 

More precisely, ground state fidelity, or simply fidelity, 
is defined here as 

F{\5) = \{X-5\\ + 5)\, (1) 

where |A) is a ground state wave function of a many-body 
Hamiltonian H{X) describing the system exposed to an 



external field A, while (5 is a parameter difference. It 
provides the most basic probe into the dramatic change 
of the wave function around the critical point and has 
the following properties 



< J"(A,(5) < 1, F{\5) = F{\,-5). 



(2) 



The recent surge in studies of fidelity follows observation 
that quantum criticality promotes its decay 0. This is 
easily justified as ground states change rapidly near the 
critical point to reflect singularities of a QPT. Therefore, 
one expects that J^(A, 5) has a minimum at the critical 
point. 

The drop in fidelity encodes not only the position of 
the critical point, but also universal information about 
the transition given by the critical exponent v. This has 
been worked out independently in the "small system" 
and in the thermodynamic [T7| limits. Broadly 
speaking the former corresponds to 5 — >■ at fixed system 
size TV, while the latter corresponds to — ?> oo at small, 
fixed 5. 

In the "small system" limit we can Taylor expand fi- 
delity in (5 0,i,[l| 



J-(A,<5)~l-<5V(A)/2, 



(3) 



where xf stands for the so-called fidelity susceptibility. 
The linear term in i5 disappears in the above expansion 
due to the normalization condition of the ground states 
whose overlap is taken (see e.g. [1]). This can be also 
seen from a linear in i5 term could make > 1 or 
simply fidelity is symmetric with respect to the 5-^—5 
transformation. Universal information can be extracted 
from fidelity susceptibility near the critical point through 
the scaling: xf(Ac) ^ iY^/'*", where d stands for system 
dimensionality. Alternatively, one may look at xf away 
from the critical point and study xf(A) ~ A^|A — Ad''''"^. 
More generally, the scaling of fidelity susceptibility is 
linked to the scaling dimension of the most relevant per- 
turbation (iGl. 
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In the thermodynamic hmit one obtains [l7j : 



A- A, 



(4) 



where A is a scahng function. We expect the above scal- 
ing resuh to hold when fidelity per lattice site (— \nT/N) 
is well defined in the thermodynamic limit, physics 
around the critical point is described by one characteris- 
tic scale of length (correlation length) given by the critical 
exponent v, and dv < 2 so that non-universal (system- 
specific) corrections to (jlj are subleading [14l - [l6j |. 
Expression (j4|) can be simplified away from the critical 
point when \S\ < |A - A^ < 1: 



lnJ-(A,(5) ~ -NS^\X~Xc\ 



du-2 



(5) 



Unlike ([3]), (H)) is nonanalytic in S near the critical point 
and the Taylor expansion of fidelity in S is inapplicable 
here. This singularity arises from anomalies associated 
with the QPT in the thermodynamic limit. All these 
results have been obtained in the zero temperature limit 
that we adopt in this work. Even though a consistent 
theory of fidelity in finite temperatures is still missing, 
several results have already been obtained [T9l - [2H . 

This article is organized as follows. Sec. describes 
theoretical approaches used to study fidelity and motiva- 
tion behind this research. In Sec. Illll we focus on fidelity 
in the XY model showing numerous analytical results. 
Sec. lyl discusses fidelity of the extended Ising model 
SHjlil. In Sec. El we illustrate some connection be- 
tween fidelity and dynamics of a nonequilibrium quench. 
Our findings are briefiy summarized in Sec. IVIl 



II. BASICS OF FIDELITY 

To start, we mention that the term fidelity is used in 
various contexts in quantum physics to describe the sim- 
ilarity between two quantum states. Therefore, it is im- 
portant to remember that we use it only to refer to an 
overlap between two ground states of the same physical 
system calculated for different values of its external pa- 
rameters. 

One of the seminal results on fidelity was obtained by 
Anderson decades ago 2j|. He showed on a particu- 
lar model that fidelity disappears in the thermodynamic 
limit. Similar behavior has been later found in other 
models and was labeled as the Anderson orthogonality 
catastrophe. More importantly, orthogonality catastro- 
phe was shown to play a role in numerous condensed 
matter systems (see e.g. [1^ and references therein). 

While being an equilibrium quantity, fidelity shows up 
in nonequilibrium dynamics of quantum systems. In par- 
ticular, one encounters it in the context of nonequilibrium 
QPTs (see [2^, [l^ for reviews on dynamics of QPTs). 
For example, the scaling of density of excited quasipar- 
ticles following a quench can be derived using it [14| - [l6| . 



Moreover, an envelope of nonequilibrium coherence oscil- 
lations in a central spin problem encodes fidelity as well 

In the context of equilibrium QPTs, fidelity has been 
proposed as an efficient theoretical probe of quantum 
criticality (the so-called fidelity approach). To appre- 
ciate its simplicity, we compare the fidelity approach to 
the traditional method of studies of quantum criticality 
based on analysis of the asymptotic decay of the two 
point correlation functions. First, both approaches re- 
quire the same input: the ground state wave functions. 
Second, they provide the same information: the loca- 
tion of the critical point and the universal critical ex- 
ponent v. Third, the conventional approach is based 
on the assumption that the two point correlations be- 
tween spins/atoms/etc. decay semi-exponentially away 
from the critical point and in a power law manner at the 
critical point [29]. The transition from semi-exponential 
to the power-law decay can be tedious to study even 
with the recent tensor network techniques (see e.g. |3C|). 
Therefore, the fidelity approach is arguably a simpler al- 
ternative at least as far as numerical calculations are con- 
cerned. 

On the analytical side, compact and accurate results 
for fidelity are scarce. Typically, models are not exactly 
solvable and so analytical results are out of reach for 
them. In exactly solvable systems situation is far from 
trivial as well. Indeed, in systems like the Ising model 
one is left with a large product of analytically known 
factors that have o nly recently been cast into a sim- 
ple expression: see and Sec. IIIII The situation is 
much different in some systems where ground states can 
be exactly expressed through finite rank Matrix Product 
States: see [3, [H, [2^ and Sec. IIVI The Hamiltonians of 
these systems, however, are "engineered" to have a pre- 
determined ground state [Q]. These difficulties motivate 
various approximate approaches. 

First, powerful numerical techniques have been de- 
ployed including tensor networks [3l| and Quantum 
Monte Carlo [lUlT^ simulations. These approaches pro- 
vide crucial insight into models that are not exactly solv- 
able (e.g. jUl), and shall be especially useful for stud- 
ies of systems with unknown order parameters, critical 
points and critical exponents. 

Second, the fidelity susceptibility approach based on 
the Taylor expansion ([3]) has been used. Simplification 
here comes from factoring out the parameter difference 
S and focusing on fidelity susceptibility that depends on 
the external field only. Despite these simplifications, it 
is typically still a mixed analytical and numerical tech- 
nique. This approach is limited to studies of fidelity be- 
tween very similar ground states whose overlap is close 
to unity. In particular, this rules out description of An- 
derson orthogonality catastrophe within this framework. 

Third, the fidelity per lattice site approach has been 
proposed [U [U [sj . It is in a sense "orthogonal" to the 
fidelity susceptibility approach as (i) it targets overlap 
between any two ground states: any S is considered; (ii) 
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it focuses on the large system limit, i.e., the opposite of 
what the fidelity susceptibility approach does; and (iii) 
fidelity typically departs significantly from unity here. 
This approach has been explored mostly numerically so 
far. 

Our studies assume two most plausible features of the 
above two approaches [l7j . First, the thermodynamic 
limit from the fidelity per lattice site approach which 
allows us to study fidelity in systems well approximat- 
ing the critical (infinite) ones. Second, small 6 as in fi- 
delity susceptibility approach allowing for derivation of 
universal scaling results such as (g]) and ([5]). Moreover, 
combination of the two assumptions allows for derivation 
of compact analytical results for fidelity in some exactly 
solvable models: the task presumably impossible other- 



Ill. XY MODEL 
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FIG. 1: (Color online). Phase diagram of the XY model ([6|. 
Critical points lie on green (thick) lines. (A) to (D) show 
different characteristic lines along which we calculate fidelity. 
The dashed "circle" separates two regions of the ferromag- 
netic phase differing by the structure of the excitation gap. 
Minimum of the energy gap is reached for fee = outside the 
"circle" (fee = vr for g < 0) and for kc = arccos[(7/(l — 7^)] 
inside the "circle" (see e.g. [ssl]). 



In this section we are going to study fidelity of the XY 
model 



N 
n=l 



1+7 



2 '^n'^n+l 



(6) 

where we assume periodic boundary conditions un+i — 
(Ti . This model is exactly solvable via the Jordan- Wigner 
transformation translating it into a free fermionic system 

0,[3ii3. 

Above g is the external magnetic field acting along the 
z direction and 7 is the anisotropy of spin-spin interac- 
tions on the xy plane. The critical exponent v for Ising- 
like critical points, g = ±1, 7 7^ and — 1 < g < 1, 7 = 
is 



(7) 



The specific case of a multicritical point located at 
5 = ±1 and 7 = needs special attention [s^. It has 
been proposed that there are two divergent characteristic 
length scales that have to be taken into account around 
it: one characterized by__the scaling exponent v = 1/2 
and the other by = 1 [37]. 
We will study here 



^ = |(3i,7i|52,72) 



(8) 



where \g, 7) is the ground state of (|6]). As depicted in Fig. 
[U we choose {gi,^i) to lie on straight lines near critical 
points/lines. 

The analytic^ expression for fidelity in the XY model 
was given in 7:]: F — Y\k>i:)\'^os(j)\/2 — 9\/i)\, where 
tan 9\, = 7i sin k/ {gi — cos k) and momenta k are given 
by pU)) . We have reworked it to the form which proves 



to be convenient for subsequent analytical studies 



k 
fk 
Pk 



k>0 

{2n + 1)t:/N, 7i = 0, 1, . . . , N/2 - 1, (10) 

(11) 



'1 1 



Pk 



= (gi - cosfc)(52 - cosfc) -I- 7172 sin^ fc, (12) 
= [72(51 ~ cosfc) - 71(32 - cosfc)] sinfc, (13) 



which is valid for any gi, 52, 7i and 72. To be more pre- 
cise, we assume even N and follow notation from [38 ] dur- 
ing diagonalization of The ground state lays then in a 
subspace with even number of quasiparticles which leads 
to ([TOl) . In that subspace the XY Hamiltonian ([6]) is diag- 

onalized to the form H = J2k>o (7fc7fc + 7lfc7-fc - l) j 
where 7^ are fermionic annihilation operators and the en- 
ergy gap is 



efe 



2 V (.9 ~ cos kY + 7^ sin^ k. 



(14) 



In the leading order fidelity © can be additionally sim- 
plified by replacing the product over momentum modes 
by an integral 



In J" = ^In/fc ~ — / dkhifk 
k>o ^ 



(15) 



which is allowed in the thermodynamic limit as all the 
integrals that we study are convergent everywhere (even 
at the critical points). We note here, however, that in 
some cases In/^ has a logarithmic singularity adding a 
subleading term to the transformation from summation 
into integration (jl5|) . This will be discussed in details in 
Sees. Hirai and lnTB] 
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A. Across the g=l critical line 

In this paragraph we foUow the path A from Fig. [TJ 
i.e., we substitute 

51,2 = l + e±S, 71 = 72 = 7, e = c\S\ (16) 

into ([8]) and furthermore assume that < 7 < 1. The 
relevant critical point is at = 1 with the critical ex- 
ponent given by 0. The relevant correlation length was 
calculated in [SJ and reads 



1 



ln|(5-\/F^T+7^)/(l-7)l 



(17) 



As long as 7^ ^ jg — 1|, it can be approximated near the 
critical = 1 by 



7 



\9 



1 



(18) 



To use this expression for the correlation length and to 
find the leading, universal behavior of fidelity we keep 
|e| <C 7^ in this section. We note that it implies that 
we stay outside of the "circle" in Fig. [TJ In this region 
the minimum of the energy gap is reached for the smallest 
value of momentum k pop . This condition also keeps us 
away from the multicritical point at 7 = and g = 1, 
which will be investigated in Sees. IIII Gl and IIIIDl 

The results of this section complement our previous 
studies from [l3| where we have considered a special case 
of the Ising model: 7 = 1. Here - besides extending our 
results to 7 7^ 1 - we discuss the problem in momentum 
space, show details of our calculations and estimate errors 
of our approximations. 

To begin the discussion, we expect the system to 
crossover from the "small system" limit to the thermo- 
dynamic limit when [l7| 



(19) 



where A denotes a point on the (g, 7) plane, while 6 is its 
displacement. The latter is in general a two-component 
vector. Above ^(A) is the correlation length in a system 
exposed to the external field A and L is the linear size of 
the system: L"^ = N. Fidelity approaches the thermody- 
namic limit result when 

mm[^{\ + 6),({\-6)]<$:L. 

On the other hand, when 

min[e(A + J),e(A-(5)] >L 

the "small system" limit is reached and the behavior of 
fidelity is dominated by finite size effects. We use the 
quotation marks to highlight that we still consider L ^ 1 
here as QPTs shall not be studied in systems made of a 
few spins/atoms/etc. 




FIG. 2: (Color online). Transition from the "small system" 
limit (left part of the plot) to the thermodynamic limit (right 
part of the plot) resulting from variation of the anisotropy 
parameter 7. Solid black line is a numerical result showing 
= \{gc — 2S, y\gc, 7)1, red dots depict the "small system" pre- 
diction (|2ip . blue crosses represent the thermodynamic result 
((27|) . The figure is prepared for 5 = 3 x 10"'^ and iV = 10®. 



We notice that the correlation length in p8)) has a 
prefactor which is linear in 7. By changing it we will 
illustrate the crossover ([T5| . To that end, we fix N, S 
and e — —\d\ {c — —1) so that one of the states is exactly 
at the critical point, i.e. J" — \{gc — 26,j\gc,^)\- The 
outcome of such a calculation is depicted in Fig. Two 
distinct regimes are visible there. In the left part of the 



plot, where 7 1, we observe In - 



-7" 



In the right 



part of the plot, where 7 ^ 1, we find InJ^ ^ ^7 • 
These two regimes correspond to the "small system" limit 
and the thermodynamic limit, respectively. 

Substituting ^ and (HH) (with fixed c) into HH) we 
obtain the condition for the crossover as 



m 

7 



1. 



(20) 



The validity of the above condition is numerically con- 
firmed in Fig. [3] This is done in the following way. 
As — In 7 is increased in Fig. [2l the slope of ln(— InJ^) 
changes smoothly from 2 (corresponding to In ~ ~7~^) 
to 1 (corresponding to In ^ —7^^). The crossover re- 
gion between the two limits is centered around 7 = 73/2 
where the local slope equals 3/2. By repeating the cal- 
culation from Fig. [2] for various system sizes N, but the 
same S, we have numerically obtained 73/2 (-^V)- A power- 
law fit described in Fig. [3] confirms that 7^/2 ^ N. Com- 
plementary analysis has been done in [l7[ where 7=1 
and (5-dependence of ([^0]) has been verified. Based on 
these two calculations we conclude that the crossover con- 
dition (PO]) holds near the Ising critical point. In fact, it 
is in an excellent agreement with our numerical simula- 
tions. 

After describing the crossover between the two regimes 
observed in Fig. [21 we can explain in detail the seal- 
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FIG. 3: Numerical study of the crossover condition (|20[) : see 
text for details. Crosses come from numerics, while straight 
line is a linear fit: ln73/2 = -15.51±0.05+(1.007±0.004) IniV 
[39I I . The fit shows that the crossover takes place when 7 ~ A*' 
as predicted by (f20)) . We used here S — 3 x 10^^ as in Fig. |2] 
and extracted 73/2 from T — Kffc — 25,'y\gc,'y)\. 



ing observed on both "ends" of the plot. We start with 
the "small system" limit. A simple analytical calculation 
based on the expansion ^ - i.e. the fidelity susceptibil- 
ity approach - shovifs that near the critical point 



J-- 1 



1672 



In ~' 7 



(21) 



because fidelity is very close to unity here. As illustrated 
in Fig. [21 this result agrees well with numerics. For com- 
pleteness, we mention that away from the critical point, 
|e| \S\, expansion ^ yields 



J-- 1 



I67R 



(22) 



in the leading order in |e| <C 7^. 

Explanation of the thermodynamic limit result re- 
quires more involved calculations. Substituting (|16p into 
(|12m3p we obtain 



Pk = (1 -I- e - cosfc)' 
Qk — 2(57 sin fc. 



7^ sin^ fc. 



(23) 
(24) 



To calculate InJ^ = \ii\{gc + e + S,j\gc + e — 6,j)\ in 
the leading universal order in S we notice that the in- 
tegral (1151) is dominated by the contribution from small 
momenta fc. This allows us to approximate (|23I24I) by 
w k'^-y'^ + e"^ — S'^ and qk ~ 2d^k. Thus, we get 




fc^7^ + 



2 V(fc^7^ + - + {25-ikf 



(25) 

Next we put (l25t into (|15l) . change the integration vari- 
able fc there to I |<5|/7, and send the new upper integra- 
tion limit, 7r7/|(5|, to 00. After all these approximations 
we obtain 



In J": 



m 

47r7 



dl In f i 



f + c^ 



1 



2Va^ 



C2 - 1 



4/2 



(26) 



The above integral can be calculated analytically, 
^""^ -\5\A{c)h, 



N 

where the scaling function 



(27) 



A(c) 



and 



1 , \c\K{c,) 

4 27r 

|c| |c|X(ci) 
4 



27r 



(|c| - l)Im£;(c2) 
47r 

_ (|c| - l)Img(c2) 

477 



for |c| < 1, 



for 



(28) 



Cl = -4|c|/(|c|~l)^ C2 = (|c| + 1)V(N-1)^ (29) 

The complete elliptic integrals of the first and second 
kind are defined as 



t/2 



\J\ — X SVC? 



r/2 



(30) 

respectively. This result is quite interesting. 

First, solution ((27)) compares well to numerics, which 
is illustrated in Fig. [5] for c— —1. Not only the relation 
InJ^ ~7~^ is reproduced, but the whole expression 
fits numerics well. We mention in passing that the same 
good agreement was obtained for 7=1 and different (5's 
in [l7| . which confirmed the general scaling prediction 
(|3]). Notice that not the anisotropy 7 but the shift in the 
magnetic field b is the relevant perturbation here. 

Second, a singularity of derivative of fidelity can be ob- 
tained from ([77)): ^lnJ^|c_^i± is logarithmically diver- 
gent (see for 7 = 1 case). To make a more transpar- 
ent connection to the former studies of Zhou and collab- 
orators who called such singularities as the pinch points 
[iSillliil^i consider -^d{gi,g2), where the scahng 
parameter is defined as 



N- 



(31) 



This can be analytically calculated from (|7f)) to be 



_d_ 

dg". 



dA 
dc 



e + S 
e/5 2-57 



27 



-A 



(32) 



where the relations between 31,2, £ and S are specified 
in (|T6)) . Above we assume gi > 52 i-e. S > 0. The 
explicit expression for (I32p is quite involved and so we do 
not show it. It predicts a logarithmic singularity when 
f/2 — >■ 1 with prefactors which depend on 51 ^ 1. More 
precisely, this singularity arises in the limit of N ^ 00, 
already assumed in ()32p . and is rounded off for finite 
systems: see Fig. |4] As will be shown below, it has an 
interesting interpretation in momentum space. 

Third, solution (|?fl) can be simplified away from the 
critical point. Straightforward expansion done for |c| 3> 1 
(but stiff |e| = \cS\ < 7^), gives A{c) ~ l/16|c| and so 



exp(-iV5Vl67|e|) 



6 




FIG. 4: (Color online). Derivative of the scaling parameter: 
-^d{gi, g2). Solid black line is the analytical result (|32[l . 
Dashed lines from top to bottom correspond to numerical 
results for = 100, 300, 1000, 3000, respectively. The field 
gi = 1.1 and 7 = 1. 



in agreement with ([5]). When the argument of the ex- 
ponent is not small, this is a new result. Otherwise, it 
coincides with expression (P^. from the fidelity suscep- 
tibility approach ([3]). 

Complementary insight into fidelity can be obtained in 
momentum space where we focus on fk ■ It is presented in 
Fig. [5] where we show prototypical behavior of fk for the 
Ising model when: (i) the ground states entering fidelity 
are calculated on the opposite sides of the critical point 
|c| < 1; (ii) one of them is obtained at the critical point 
|c| = 1; (iii) they are both obtained on the same side of 
the critical point |c| > 1. 

As illustrated in Fig. [5l if Afc = 2t:/N > \S\/^ only 
fk ~ ^ contribute to fidelity, and we end up in the 
"small system" limit. The system is too small to monitor 
changes of fk between zero and unity. StilL we are able 
to observe universal finite size effects |ll| - [l6| in (I^Tt . 
In the opposite limit of Afc ^ |<5|/7, we are able to ap- 
proach the exact (TV — )■ 00) thermodynamic limit and 
the product ^ is well approximated by the integral ([15]) 
[but see the discussion below as well]. Now momenta k 
are dense enough to monitor leading changes in fk be- 
tween zero and unity. 

Still, to see sharp nonanalyticity (i.e. the pinch point) 
in fidelity when one of the states passes through the crit- 
ical point we need to have iV — >■ c» as demonstrated in 
Fig. m It is caused by discontinuity of fk at zero mo- 
mentum when |c| = 1. Indeed, from (|lll23l24p it is easy 
to calculate the limits 



limfe_ 
limfe_ 

limfe- 



■k^ fk = 0^ 

■fe. 



Jk 
fk = l 



for |c| < 1, 
for |c| = 1, 
for IcI > 1, 



(33) 
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FIG. 5: (Color online). Thermodynamic vs. "small system" 
limit near the Ising critical point. Solid lines correspond to 
fk calculated for 7 = 1, 5 = 0.01, and - from bottom to 
top - c equals 0, 1 and 1.5, respectively. Small red crosses 
mark fkS contributing to fidelity calculated for A'^ = 3000. 
Large black crosses show the same, but when the system size 
is A'^ = 50. For = 3000 the thermodynamic limit is reached. 
A' = 50 represents the "small system" limit. Notice that in 
the "small system" limit the resolution in momentum space is 
insufficient to account for abrupt changes of fk taking place 
at small momenta. 



happening on the momentum scale given by the inverse 
of the larger correlation length scaling as \5{\c\ — l)|/7, 
which is divergent at |c| — >■ 1 (notice that we have two 
correlation lengths here for the two ground states enter- 
ing fidelity). It results in rounding off of the derivative of 
the scaling parameter d — — hiF/N around the critical 
point for finite N due to insufficient sampling of fk near 
kc — Q- We also notice that the approximation (|25|) that 
we make to obtain (|27p correctly capture this discontinu- 
ity. 

Finally, we can discuss errors resulting from our ap- 
proximations. On the one hand, approximating exact 
expression for fk (|lll23l24p with ^ in the integral ^ 
leads to 0{5'^ /^^) error in (|27p. This calculation is quite 
technical and has been deferred to the Appendix. On the 
other hand, there are errors connected with estimation of 
the product ^ by the integral Main contribution 

here comes from the logarithmic divergence of In fk at 
k = Q when 

\c\ < 1. 

It results in a subleading shift which has to be added to 
the right-hand side of (fT5|) . This shift saturates to In 2/2 
when the size of the system TV is large enough to well 
sample the logarithmic singularity. It happens when the 
system size iV is much larger than the larger of the two 
correlation lengths proportional to 7/|(l — \c\)5\. It can 
be seen, e.g., when we expand fk to the lowest order in 
k around kc = 0, 



where kc = here. This is presented in Figs. [S] and H) 
In the latter, it is shown that the jump of fk is roughly 



fk 



jk 



(34) 
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FIG. 6: (Color online). Discontinuous change of fk around 
c = 1. The curves from bottom to top correspond to c equal to 
0.98, 1 and 1.02, respectively. Red crosses, similarly as in Fig. 
m mark fk's contributing to fidelity calculated for N — 3000. 
The parameter shift 5 — 0.01. Lowest order Taylor expansions 
of fk are listed (but not plotted) for all three curves. 



for |c| < 1. Using this result we obtain 



0<fc<fccutoft 



fecutott 



dkln{ak) ___ 



Above, a — 7/|5|(l — c^), fc's are given by dTO 



(35) 
and 

^cutoff restricts summation/integration to small momenta 
for which the expansion is meaningful. We get the 
equality in (1351) in the limit of fccutoff / Afc — >■ oo when the 
logarithmic singularity is sampled densely. This can be 
easily obtained using the Stirling formula. 

Including the above discussed errors/corrections into 
(l27t we obtain 



c| < 1, 
c| > 1, 
(36) 



lnj;_/ -\S\A{c)/j + \n2/2N + 0{S^/-f^) for 
N " \ -|<5|A(c)/7 + ©((52/73) for 

which is equivalent to 



V2exp [-N\S\A(c)/j + NO{S^/-f^)] 
exp [-N\5\A{c)/-f + NOiS^f-f"^)] 



for |c| < 1, 
for |c| > 1. 
(37) 

Note that the difference between |c| < 1 and |c| > 1 
cases comes from the lack of a logarithmic singularity in 
In fk in the latter. As illustrated in Fig. [71 numerics 
compares very well to (|36p and (|37)) . The prefactor \/2 is 
seen in numerical simulations when iV 3> 7/|(l — |c|)5|. In 
derivation of and ([57]) we have neglected other errors 
coming from changing the product ^ into the integral 
([15]). Those corrections are present around c — ±1 but 
are subleading and disappear in the limit of oo (e.g. 
they smooth out the pinch point nonanalyticity) . 

Finally, we would like to stress that a subleading cor- 
rection to the argument of the exponent in (j36p . the 
In 2/2 A^ term, increases fidelity by the ^2 factor ([571) . 
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FIG. 7: (Color online). Panels in the left and right column 
illustrate In and InT /N\S\ as a function of N, respectively. 
From top to bottom rows show results for c equal to 0.9, 1 and 
1.1. Black solid lines show numerics. Red dashed lines present 
— ^(c)A|5| in the left column and —A{c) in right column. 
Green dashed-dotted line in the top row shows — ^(c)A'^|(5| + 
In 2/2 in the left panel and ~A{c) + ln2/2iV|5| in the right 
one. See discussion around p6p and (|37p for more information 
about the logarithmic correction. The parameter shift S = 
10"^ and 7=1. 



Its influence can be neglected when we study the scaling 
parameter, —\nT/N from (15^1) . instead of fidelity. This 
is presented in the right column of Fig. [71 



B. Across 7=0 critical line 

In this paragraph we follow the path B from Fig. [1] 
i.e., we substitute into ^ 



71,2 = e ± (5, 31,2 = 5 e (-1, 1): e = c\S\. 



(38) 



Two remarks are in order now. First, we shift 7 instead 
of g here, but use the same symbols, S and e, for the 
shift as in Sec. IIII Al Therefore, some care is required to 
avoid confusion. Second, critical exponents z and v are 
the same as for the Ising universality class studied in Sec. 
IIII Al Other critical exponents, however, differ between 
the two transitions j35| . 

Expanding (|17p in small 7 relevant for our calculations, 
the correlation length reads 



when 1 ^ 7 < a/1 — (jr2_ This can be used to predict the 
crossover condition from ([TOl) . We focus here on large 
enough N to study the thermodynamic limit. Further- 
more, we restrict our calculations to |e|, |<5| ^ \/^ ~ 9^ 
in order to make the problem analytically tractable. 

The crucial difference with respect to the calculations 
from Sec. IIII Al is that at the critical points - lying on 
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FIG. 8; (Color online). Plot of fk for the transition across 
the 7 = line We choose g = 0.5 and 5 = 0.1 here. 

The curves from bottom to top correspond to c equal to 0, 1, 
and 1.5, respectively. Dominating contribution to fidelity is 
centered around kc = arccos(gi) — n/S. 



the 7 = and g e (— 1, 1) line - the gap in the excitation 
spectrum (jl4p closes for momentum 



kc = arccos((7). 



(39) 



As can be expected, fidelity in the thermodynamic limit is 
dominated by contributions coming from momenta cen- 
tered around kc (Fig. IH). 

Putting ([38|) into (|12|13|) one gets the following exact 
expressions 



Pk 



[g - cos kf + (e^ - 5^) sin^ fc, 



qk = 2d{g — cos k) sink. 



(40) 
(41) 



We expand the above around dominant kc to get pk ~ 
{l-g^){{k-kcy+e^-6^) and ^ (l-g^)2^( fc-fc,) , 

where we have utilized the condition |e| <C \/l — 5^. 
Thus, we get 



In/fc 



im i 
2 1 2 



(fc-fee)2+. 



2\/((fc - fce)2 + £2 _ 52)2 + 4^2^^ _ fc^)5 



(42) 



We integrate it over k ((TS]), parameterize e as c|i5| (note 
that e is now the shift of the anisotropy parameter), 
change the variable k in ([T5|) to fee + Z|<5|, and send the 
new integration limits, —kc/\6\ and (tt — fcc)/|<5|, to =Foo 
respectively. 

After all these approximations we arrive at 



In J"sr 



m 

4-K 



dl In \ - 



+ - 1 



2^(Z2+c2-l)2+4P 



(43) 

We have introduced above the symbol J^3mooth because 
it turns out that for |c| < 1, i.e. when the two states 
used to calculate fidelity are obtained on opposite sides 
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FIG. 9: (Color online). Comparison between the smooth part 
of fidelity (|44p - given by the red dashed line - and the exact 
numerical result (solid black line). In the large A'^-limit and 
for |c| > 1 oscillations of fidelity are absent, which is presented 
here for c = 2 (top panel). At c = 1 some oscillations of fi- 
delity are visible for not-too-large-systems (middle panel), but 
they gradually disappear for large A^. For |c| < 1 pronounced 
oscillations of fidelity surviving in the thermodynamic exist: 
bottom panel prepared for c = 1/2. The plots are made for 
5 = 0.002 and g = 0.99. 



of the critical point, there will be oscillatory corrections 
to fidelity (Fig. [9|). Before describing them, we analyze 

smooth- 

Except for the integration range, the above integral is 
exactly the same as the one for the Ising model in (|26p. 
Therefore, we obtain 



In J's, 



DOth 



N 



-\5\2A{c) + O {5^ 



(44) 



where A{c) is defined in ([^5)) . We estimated numerically 
the error between the exact integral (jisp and the approxi- 
mated one (|43p . what is presented in the Appendix. Note 
that (HI]) quantifies decay rate of fidelity with system size 
even when the oscillations are present (Fig. |9]). 

As pointed out in e.g. (35| . the anisotropic transition 
has the same critical behavior as a pair of decoupled Ising 
chains. This can be seen also through fidelity: there is a 
prefactor of two in (|44p absent in ((27)) from the previous 
section. Looking from the momentum space perspective, 
the factor of two conies from different location of the 
momentum kc- For the Ising model considered in Sec. 
nil Al kc — and only k > kc contribute to fidelity ([9]). 
Here fcc 7^ and so both k > kc and k < kc add up to 
fidelity doubling the result [additivity of this effect can 
be seen from p5|) ]. 

We mention also that the pinch point singularity is 
predicted by (|44p . It is again visible in momentum space 
through discontinuity of fk and described by p3p : notice 
that now fc^ 7^ This is iUustrated in Fig. [S] 
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FIG. 10: Plot of /fe (fTTj) for the transition across the 7 = 
line for S — 0.01, c = and several values of g such that (a) 
fcc = ^3§§§, (b) kc = 7r/4, and (c) kc = tt^^ [note that kc 
depends on g (|39[) ] . Pluses mark discretization in momentum 
space for the size of the system N = 2000. We see a systematic 
shift of /fc's contributing to fidelity as we slightly change g. 
This results in oscillations of fidelity from the lowest panel of 
Fig. [51 Parameter from (|45|l equals 1, and 1/2 for panels 
(a), (b), and (c), respectively. 
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FIG. 11: Schematic plot illustrating the shift (p from (I45|) . 
Crosses show discretized momenta, the dot denotes position 
of kc (|39|) . the dashed fine goes through the mid point between 
two discretized momenta being closest to fee = (2m + <j))-n /N . 
Note that > (0 < 0) when kc is to the right (left) of the 
dashed line. Naturally, when (j> approaches ±1, kc collapses 
onto one of the discretized fc's. 



k mi- 



Nkr 



modulo 2 e (-1, 1] 



(45) 



Its definition is illustrated in detail in Fig. [TTJ The ex- 
amples oi (j) — 1, and ^ are presented in Fig. [TOl To 
extract the oscillatory part we focus on contribution to 
fidelity coming from momenta around kc- Taking the 
lowest order term of the Taylor expansion of (lll[) around 
kr we obtain 



fk 



\k kr 



(46) 



Now we are ready to describe the oscillatory correc- 
tions to (j44p depicted in Fig. |9l These effects, visible 
for 

|c| < 1, 

are washed away by the above continuous approximations 
relying on the limit of — >■ c» that overlooks discretiza- 
tion of momenta k. The reason behind the oscillatory 
behavior is easily identified by looking at Fig. [TU] In- 
deed, fidelity is a product of /fc's ^ and as such it is 
sensitive to the fk factors approaching zero. For |c| < 1, 
/fc's stay close to zero when k k kc (as shown in Fig. 
[5] this is not the case for |c| > 1). The arrangement of 
fc's around kc - fixed by discretization (fTO)) - determines 
fidelity oscillations. 

Qualitatively, fidelity reaches its minimum at zero 
when one of fc's equals exactly kc- This happens when 
kc = (2n + l)7r/iV for some integer n £ [0, 1]. Then 

fk=k^ ~ and a single momentum mode takes fidelity 
down to zero even for finite N. This case is illustrated in 
Fig. [TOk . On the other hand, when kc lays symmetrically 
between two discretized fc's, Fig. [TUb . we can expect that 
fidelity is near its maximum. All other "orientations" of 
discretized momenta fc with respect to fcc, including the 
one depicted in Fig. [TUb . add to the oscillation presented 
in the bottom panel of Fig. [3] 

Quantitatively, we introduce the parameter (j) measuring 
distance of fcc from the mid point between two momenta 



where a — 1/\S\{1 — c^) is introduced for convenience 
[compare with (|34l)]. Without going into details, we de- 
fine fccutoff estimating the largest |fc — fcd for which (j46p 
provides a reasonable approximation to the exact expres- 
sion for /fc. Following notation introduced in Fig. [TTl we 
easily find that fk coming from the m-th momenta to the 
left (plus sign) or right (minus sign) of fcc equals 

/fc«a^(2TO-l±(/.). 

Now we can find how the asymmetry of momenta fc 
around fcc modifies fidelity. We do it by factoring out 
the contributions coming from the (f) = case where no 
asymmetry is present. Namely, we calculate 

Ilfc/fc _ " TT°" [-, _ ^ 
UkhU=o-J,{ [ (2m- 1)2 J 2- 

(47) 

The last equality is reached when fccutoff /tt ^ 1- The 
change of the product into the cosine can be found in 

To use (|47|) . we notice that for = fidelity is 
given by 2jsinooth (|44p . Thus, for any (j> it is given by 
2 cos(7r(/)/2) ^"smooth as long as |c| < 1. The origin of the 
prefactor of 2 which multiplies J^3lnooth is the same as in 
Sec. nil Al We have to correct lnJ^5mooth by a sublead- 
ing term of In 2 to account for logarithmic singularities of 
In/fc near fcc ((35|) . Note that for = the displacement 
of momenta fc around fcc is equivalent to displacement 
of momenta fc around fcc = for the Ising model - Sec. 
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FIG. 12: (Color online). Plot of fidelity for the transition 
across the 7 = hne for 5 = 0.002 and various parameters 
g and |cj < 1. Black line and crosses show numerical results. 
Red dashed line and circles show theoretical prediction (|48l) . 



IIII Al Thus, we can write our final prediction for fidelity 
in the following form 



2 I cos 
exp f- 



(^)jexp [~N\5\2A{c) 
-N\S\2A{c) + NO{5^)] 



NO{S^)] for |c| < 1, 
for |c| > 1, 
(48) 

where was employed to simplify the oscillating fac- 
tor. The factor of two and sinusoidal shape of oscilla- 
tions is reproduced by numerics when the system size 
A'' is much larger than the larger of the two correlation 
lengths proportional to 1/|(J(1 — |c|)| (compare to Sec. 

MM- 

The oscillating part of p8)) is illustrated in Fig. [121 As 
we see in panels (a) and (b), an agreement between (j48|) 
and numerics is remarkable. The limits of applicability of 
(|3H|) are illustrated in panel (c) of Fig. [T2] The left part of 
this panel shows discrepancies resulting from application 
of pS)) to systems whose size N is not much larger than 
the larger correlation length (~ 1/|(5(1 — |c|)|) -still notice 
that the period of oscillations is predicted correctly. The 
right panel reveals that for large enough system sizes, 
and the very same parameters, the perfect agreement is 
recovered. We note in passing that those oscillations will 
survive in the "small system" limit making the fidelity 
susceptibility approach (jS]), which assumes that fidelity 
stays close to unity, inapplicable to this problem. 

We expect that the same mechanism as described 
above should be able to explain oscillations of fidelity 
(peaks in fidelity susceptib ility ) observed in the gapless 
phase of the Kitaev model [41[ . 

Finally, it turns out that we can extend the results 
of this section. The condition |<5|, |e| <C \/l — g^, which 
we have been using until now and which is necessary to 
justify approximations used above, turns out to be too 
strong. Numerics shows that equation pS)) holds for any 
\S\, |e| ^ 1 and —l<g<l when N is large enough. 



0.995 



71/4 



jr/2 



37t/4 



FIG. 13: fk along the Ising critical line (|49|) for e = 1 and 
S = 0.05,0.1,0.2 (top to bottom). The minimum of ft is 

reached for k = arccos (^j^^t^jj^ — arccos ^ \^1'a ^ • 



C. On the g=l critical line 

In this paragraph we follow the path C from Fig. [1] 
i.e., we substitute into (|S]) 



7l,2 



91.2 



1. 



(49) 



This choice of the path is much different from what we 
have considered until now. Indeed, both ground states 
are calculated at points on the XY phase diagram that 
lie on the critical line. In particular, it means that the 
correlation length at any of the points is of the order of 
the system size N. 

We start by discussing fidelity calculated away from 
the multicritical point {g — — 0). This requires 

\S\ « \e\ < 1. 

Notice that e is not assumed to be small here. 

The exact expression for fidelity is given by substituting 

(|49| into ([MTS]). It can be simplified as follows. Since 



Pk (1 - coskf + (e 
Qk = 25(1 — cos A:) sin fc. 



2 S^) sin^ k, 



(50) 
(51) 



for the C path, we see that pk is positive and q'^/p'^ < 
S'^ /{e^ — S^) is small. Expanding In fk in qt/pk we obtain 



It can be analytically integrated over k from to vr (|15l) , 
and then expanded in 5 to yield 



In J" 



8kl(i + kl)2 



(52) 



Equation ()52|) provides an interesting result showing that 
the universal part of fidelity, present in Sees. IIII Al and 
IIIIBI is absent here. Therefore, only the nonuniver- 
sal part proportional to quadratic displacement of the 
ground states is present. 

Looking at fidelity from the momentum angle. Fig. 1131 
we notice that /^'s are very different here from what we 
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have studied in the previous sections. Not only f^s stay 
close to unity for all momenta, but also the main con- 
tribution to fidelity comes in general from large fc's that 
are not associated with disappearance of the gap: the 
gap (|T4l) closes at fcc = on this line. This shows that 
the ground states used to calculate fidelity differ only by 
the nonuniversal contribution away from the multicritical 
point. This is to be expected from the renormalization 
group perspective because we calculate fidelity between 
states which differ only by an irrelevant variable. 

Now we would like to show what happens near the 
multicritical point and so we assume that 

|eU^I«i, 

and use the familiar parameterization e = c\S\. The dom- 
inant contribution to fidelity comes now from momenta 
around kc — 0. 

Approximating pfc « fc^(fc^/4-|-e^ — (5^) and qk ~ k^S from 
(|50|51|) and using the same methods as in Sec. IIII A| we 
obtain 



lnJ" = 



- N[\S\2A{c)+0{6^)] 

- N [\S\2A{c) + 0(6^)] 



In 2/2 for |c| < 1, 
for |c| > 1, 
(53) 

where the scaling function A{c) is given by (|28|) . 
This is a pretty interesting result showing that despite the 
fact that the ground states entering fidelity are calculated 
on the critical line, there is a universal contribution to 
fidelity. This contribution can be traced back to the mul- 
ticritical point: the dominant critical point on the g = I 
line. Evidence for that is twofold. First, the universal re- 
sult (j53p appears as we approach the multicritical point. 
Second, the scaling exponent v, that can be extracted 
from (f53| by comparing it to ([4]), equals 1. It matches 
the predicted value for this multicritical point along the 
g = 1 line/direction [37|- It should be also mentioned 
that the concept of the dominant critical point has been 
introduced in |42l] in the context of quench dynamics. 

Finally, it is interesting to see how ([55)) crosses over 
into ([52)) as we go away from the multicritical point (the 
limit of |c| S> 1). Assuming that \S\ <^ |e| <^ 1, equation 
([5^ can be simplified to 

In J- 6^ 



N 



which matches ([5^ in the leading order in e. Note, 
however, that one cannot apply ()53p to arbitrary dis- 
tances away from the critical point because it is valid for 
|e| = |c<5| <^ 1 only. For larger e's non-universal contri- 
butions dominate fidelity ([5^ . 



D. Near 7=0 and g=l: multicritical point 

In this paragraph we follow the path D from Fig. [TJ 
i.e., we substitute into ((5)) 

51,2 = l + e±S, 71,2 = a(e ±6), e = c\6l (54) 



-15 



-20 
^ -25 



-30 


^ -10 

Oh 

■S -20 
"a 

^ -30 
-40 



- (a) 


1 ' 1 ^ 


In F ~ -N^ 




'-\^ 


■""^^ In F ~ -N 


6 


InN 12 



■ (b) 




In F ~ -5^ 










1 , 1 



-21 



In |8| 



-14 



FIG. 14: (Color online). Transition from the "small system" 
limit to the thermodynamic limit near the multicritical point 
located at g — 1 and 7 = 0. (a): the transition as a function 
of the system size A'^ at fixed S — 10~*. (b): the transition 
as a function of parameter difference S at fixed A*" — 3000. 
In both plots solid (black) line shows fidelity for c = 1, i.e., 
T = 1(1, 0|1 + 25,2(5)1, while the dashed (red) line shows the 
c = 5 case, i.e., T = \{l + i5,4:S\l + 65,6S)\ (O. Note that 
a = 1 is assumed here. 



where < a < 00 is the slope of our path. We restrict our 
studies to gi,2 > 1 and 71^2 ^ 0: both ground states used 
to calculate fidelity are obtained on the paramagnetic 
side (Fig. [1]). This requires e > \5\ or equivalently c > 1. 
Fidelity calculated along this path is "influenced" by the 
multicritical point at gc — 1 and 7c = 0. 

The condition for the crossover between the "small sys- 
tem" limit and the thermodynamic limit near the multi- 
critical point reads ([T9)) 



ivi^r - 1. 



(55) 



As we deal here with a multicritical point that is charac- 
terized by more than one divergent length scale, it is im- 
portant to carefully verify the above prediction to make 
sure that the relevant v is used in (I55p . 

The crossover is illustrated in Fig. [M) where distinct 
scalings of fidelity with either the system size or the 
parameter difference 6 are easily observed. In Fig. [Hb 
the parameter difference S is kept constant and the sys- 
tem size N is varied. We see that for small system sizes 
InJ^ ~ —N"^, while for large ones InJ^ ^ ~N. In Fig. 
[Tib the system size N is kept constant, while the pa- 
rameter shift 5 is varied. Again two regimes appear: for 
small 5 we have InJ^ ~ — <5^, while for larger 5 we get 
lnJ-~-|(5|3/2. 

The location of the crossover can be studied in exactly 
the same way as in Sec. IIII Al Briefly, as the system size 
N is increased in Fig. [T4k . the slope of ln(— In J^) changes 
smoothly from 2 (corresponding to InJ^ ^ —N"^) to 1 
(corresponding to In ~ —N). The crossover region be- 
tween the two limits is centered around N — where 
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FIG. 15: (Color online). Numerical study of the crossover 
condition (|55p . Crosses come from numerics (see text for de- 
tails), while straight lines are linear fits Is^. In both panels 
upper (lower) data sets correspond to c = 1 (c = 5). Panel 
(a): we fit a and 6 coefficients of In = a + 61n5. For c = 1 
(c = 5) we obtain h = -0.499 ±0.0003 (6 = -0.4992 ±0.0005). 
Panel (b): we fit a and h coefficients of In (57/4 = a + b\nN. 
For c = 1 (c = 5) we obtain h = -1.989 ± 0.005 (b = 
— 1.9987 ± 0.0008). The numerics for these plots is done with 
a = 1. The fitting coefficient a is not listed as it is of minor 
interest. 



the local slope equals 3/2. By repeating the calculation 
from Fig. [T3]for various (5's we have numerically obtained 
-^3/2 (i^)- A power-law fit described in Fig. [TSb reveals 
that supporting v ~ 1/2 in (1551) . Sim- 

ilar analysis can be performed on data from Fig. 114b . 
Indeed, now we look at (57/4 (TV) such that the slope of 
ln(— In J^) equals 7/4, i.e., is half-way between 2 and 3/2 
observed on both "ends" of Fig. [Hb . Fig. [TSb . where 
the power-law fit is performed, shows that (^7/4 ~ iV~^ 
again pointing to = 1/2. Therefore, we conclude that 
the crossover condition ([55l) holds near the multicritical 
point with v = 1/2 on the paramagnetic side. 

We focus on the thermodynamic limit again. Substi- 
tuting (inil) into (|12I13P we obtain 



Pk = {l + e-cosky -S^ 
Qk — 2ch5(1 — cosfc) sinfc. 



2/ 2 

a e 



(5^) sin^ fc, 



When both points - (51, 71) and (32, 72) - are on the 
paramagnetic side, then pk > and ql/pj. < ci\S\/2 ^ 1 
is a small parameter in which we expand In fk ■ This 
brings us again to 



In/fc 



(56) 



Integrating ((56|) over k from to vr p5|) . and then ex- 
panding the resulting expression in S we get 



In J" 
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FIG. 16: (Color online). Comparison between numerics (solid 
lines) and analytics [pluses given by (|57|) ]. In all panels A*' — 
10^ and a — 1. The parameter shift 5 is chosen to be small 
enough to keep the C'(|(5|^''^) corrections to H57p negligible. 



The scaling function for the multicritical point is given 
by 



A 



MCP 



(c) 



(c ± 1)3/2(3 - 2c) + (c - 1)3/2(3 ± 2c) 
16V2 



-|,5|3/2a2^MCp(c) + S^a^A + 0{\d\^^^). (57) 



which is valid for c > 1. This result is illustrated in Fig. 
1161 where small enough 6 has been chosen to keep the 
terms OdfJj'^/^) negligible. We assume such choice of S 
from now on. The result (1571) is interesting for several 
reasons. 

First, the exponent 3/2 does not fit the expected diy 
value from (|4]): neither i/ = 1/2 nor v = 1, both con- 
sidered near the multicritical point (37j . agree here. Ex- 
planation of that anomaly is beyond the scope of this 
work. 

Second, due to remarkable simplicity of (j57|) we can an- 
alyze in detail the interplay between the nonanalytic and 
the subleading contribution to fidelity. The nonanalytic 
term dominates over the subleading one when 

.Mcp(c) » 5'aV4 \S\ « 16Amcp(c)' < 1/4. 

Additionally, the "thermodynamic limit" condition 
N\/S 1 must be satisfied so that ([57]) holds. If these 
two conditions are fullfield. In J" « —N\S\^/'^a'^AMCp{c)- 
This does not, however, guarantee that T is well approx- 
imated by exp(— A^|(5|3/2Q;2^j^gp((;)) Indeed, the latter 
requirement puts an additional bound on the parameter 
shift S: 

< 1 =^ exp{NS^a'^/4) w 1. 

All these conditions can be simultaneously satisfied. 

Third, equation (|57|) shows that fidelity near the mul- 
ticritical point - at least in the paramagnetic phase dis- 
cussed here - does not have to be small in the thermody- 
namic limit. This is in stark contrast to what we found 
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IV. EXTENDED ISING MODEL 
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FIG. 17: (Color online), 
space angle that In ~ 



We show here from the momentum 
-|(5|^/2 dSZl). We plot the rescaled fk 
functions for c — 1 and a — 1: one for S = 0.001 (solid black 
line) and the other for 5 — 0.01 (dashed blue line). As the two 
curves overlap very well, we see that (I — fk)/\S\ — h{k/ ^/\5\), 
where h is some non-singular function. This implies that 
InJ-" ~ J dkln fk ^ — [fk ~ 1 so we approximate 

In/fe ~ fk — 1). This confirms leading scaling of (|57p with S 
and independently confirms the condition (|55|) with = 1/2. 
We also note that these rescalings work in the limit of 5 — >■ 0. 



near the Ising critical point, where the thermodynamic- 
limit condition (pHl) implied smallness of (l?f)) . Here, wc 

can have A^\/|5[» I and stiW eyip{-N\S\^/^a'^ Amcp{c) + 
iV(5^a^/4) anywhere between and 1. 

Fourth, similarly as in the Ising chain one can easily 
derive from (1571) transition to the analytic limit in 6 far 
away from the critical point. Considering e ^ \S\ or sim- 
ply O I (but still e = c\6\ < 1), Amcp(c) =i 5/32 V2c 
and so 



J" ~ exp -J^TV 



5a' 



32V2e 



(58) 



This is again a new result. If we expand (j58p assuming 
that the argument under the exponent is small, we get 
a fidelity susceptibility-like expression in the thermody- 
namic limit 



5a' 
16%/2^' 



In this section we consider the Hamiltonian Q 

N 

^ = E [-2(1 - 9')«+i - (1 + 9)'< 



n=l 



+ (1 - 5)'«+l<4 



(59) 



with periodic boundary conditions Un+N = cr„. Its 
ground state is given exactly by finite rank Matrix Prod- 
uct State [9']. It has the critical point at gc = 0, which is 
characterized by the critical exponents v = 1 and z = 2. 
In contrast to the Ising model, neither entropy of en- 
tanglement nor the ground state energy is singular at 
this critical point. Still, as we will see below, the basic 
features of that transition are captured by fidelity. Intro- 
ducing the notation \g) for the ground state of ([5^ . and 
defining J- as |(5i|52)| one obtains that 



(1 + + (1 - Vam) 



N\ 



+ (1 - 9ir^{i+92r + {1-92^ 

This is an exact expression from [2^ [23| . We are going 
to simplify this result to discuss it from the perspective 
relevant for our approach. 

We follow notation from Sec. IIIII writing 



.91 



+ 92 = f - S, e = c\S\. 



(61) 



We stay close to the critical point so that \S\, \e\ ^ 1 and 
approximate (|60p with hyperbolic functions 



J": 



I cosh {NVe^~S')\ 



^cosh {N{S + e)) cosh {N{e - 6)) 



(62) 



This simple expression allows us to easily study the 
"small system" limit and the thermodynamic limit. 

First, we focus on the "small system" limit, say (5 — >■ 
at fixed N for simplicity. Expansion of in S, the 
fidelity susceptibility approach (jSj, results in 



J-- 1 



1 



cosh^(eA^) 



+ 



tanh(eA^) 

7n 



(63) 



This expression depends on two parameters: N\S\ and 
iV|e|. As we keep the former small, we can focus on its 
7V|e| dependence. When |e|iV ^ 1 then ([63)) simplifies to 

J" ~ 1 - S^N' 



The above equation can be also obtained when we start 
with ([3]) and expand the resulting expression for xf in 
e. It predicts that fidelity susceptibility away from the 
multicritical point in the paramagnetic phase scales as 

Finally, we mention that one can find the |(5p/^ scaling 
of InJ-" by looking directly at fk- This is explained in 
Fig. [HI 



reproducing expected scaling of fidelity susceptibility at 
the critical point 11-16]. In the opposite limit of |e|iV ^ 
1, (1551) can be written as 

T c^l-d'N/2\e\. (64) 

This again agrees with the scaling predictions [TT| - [l6| . 
Finally, note that when (p5|) holds, fidelity stays very 
close to unity. 
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Second, we concentrate on the thermodynamic hmit, 
say iV — >■ oo at fixed S, where we obtain from (p^ 



'1-C2 

J" - < \/2cxp(-|(5|iV) 
exp(^\S\N 



cxp(-|(5|iV) for |c| < 1, 
for |c| = 1, 
- l)^ for |c| > 1. 

(65) 

To be more precise, we mention that (|65p pro- 
vides a good approximation to the exact result when 
max [N\e + 5\, N\e — 5\] ^ 1, which is a thermodynamic 
hmit condition based on (fTO]). 

Similarly as for the XY model, we introduce the scaling 
function Amps{c) as 



A 



MPS 



1 



for |c| < 1, 



|c|- 



1 for Id > 1. 



(66) 



Amps{c) is nonanalytic when |c| = 1, i.e., when one of 
the states is exactly at the critical point, signaling the 
pinch point singularity [H, [3l|, [s^]. As expected, this 
singularity is rounded off in finite systems. 
When |c| :s> 1, or equivalently \6\ <^ |e| <^ 1, then the 
scaling function Amps approaches zero as l/2|c| leaving 
us with 

T c^exp{-S^N/2\e\) . 

This coincides with the fidelity susceptibility result (|64|) 
when the argument of the exponent is small, but provides 
a different result otherwise. 

Now we will study the origin of oscillations and pref- 
actors appearing in (|65|) illustrating how the analytical 
techniques from Sec. IIIII can be applied to other mod- 
els as well. Hamiltonian ([59]) can be diagonalized us- 
ing exactly the same formalism as in Sec. IIIII Indeed, 
the Jordan- Wigner transformation translates ([59l) into a 
chain of noninteracting fermions which can be solved us- 
ing Fourier and Bogoliubov transformations. We assume 
even N and follow notation from during diagonal- 
ization of (f59| . The ground state lays then in a sub- 
space with even number of quasiparticles. In that sub- 
space the Hamiltonian (l59l) is diagonalized to the form 

H = Efe>o ffe {ihk + llk^~k - l) , where are the 
fermionic annihilation operators and the energy gap is 
given by 



ek=4{l + g^-{l-g^)cosk) 



for kc = 0. For 
8q2 + 2P, which 



It reaches zero at the critical point gc 
small k and g we can approximate 
confirms the critical exponents v — 1 and z = 2 (at the 
critical point the gap scales as fc^, while near the critical 
point it closes as \g — gd^")- 

Fidelity is given by J" = nfc>o I cos(6'^/2 - 6*2/2)1, 
where this time 0^'s are defined by 



tan( 



2{1~ gf) sink - {1 - g,)"^ sm2k 
(1 + g,Y - 2(1 - 5,2) cos fc -I- (1 - g,Y cos 2k ' 




FIG. 18: (Color online). The lines correspond to fk plotted 
for 5 = 0.01 (|69|) . From bottom to top c equals 0, 1 and 1.5, 
respectively. For c = the ground states entering fidelity are 
calculated on the opposite sides of the critical point, for c = 1 
one of them is obtained at the critical point, while for c = 1.5 
they are both obtained on the same side of the critical point. 



After some algebra we arrive at 

^ = X{\h\, 



(67) 



fc>0 



k (2n+l)7r/iV, n = 0, l,...,iV/2-l, (68) 
Vk 



fk i—ry 

Pfc = 1 + .91.92 - (1 - 5132) cos fc, 
(Ik = (.91 - .92)sinfc. 



(69) 

(70) 
(71) 



For the sake of presentation we allow fk to take nega- 
tive values here and use the parameterization (|61|) . The 
prototypical behavior of fk is shown in Fig. [181 

We expand pk and qk around kc — and for |e|, |^| ^ 1 
we get Pk ~ fc^/2 + 2(e2 - ^ 2Sk. Thus we 

can approximate 



In I A 



1 



In 



(fc2/4- 



(72) 



Before we discuss fk in details, we start the calcula- 
tion of fidelity by approximating the product in (j67p 
by an integral. Similarly as in Sec. IIII Bl we introduce 
InJ^sinooth = ^ dk\n\fk\. In order to calculate it we 
use approximation ()72|) . change the variable k into \S\l 
and send the new upper integration limit, 7r/|(5|, to infin- 
ity. This yields 



In J", 



smooth 



N r 



dl In 



(/2/4+(c2-l))2 



/O vGV4+(c2-l))2+/2 

The last integral can be calculated analytically: 

In J' gmooth 



N 



-\S\AMPs{c) + OiS'), 



(73) 



where Amps{c) is defined by ([Mjl . We notice that 
smooth (continuous) representation of fidelity reproduces 
leading exponential decay from (|65l) . It also agrees with 
the general scaling results (|H) and ([S]). The 0{S'^) error 
in the integral can be verified numerically in a similar 
way as described in the Appendix. 
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FIG. 19: (Color online). Discontinuous change of around 
c= \. The curves from bottom to top correspond to c equal to 
0.98, 1 and 1.02, respectively. The parameter shift & = 0.01. 
Lowest order Taylor expansions of around fee = are listed 
(but not plotted) for all three curves. When |c| < 1 then 
crosses zero for fco / 0. 



As noticed earlier, Amps{c) is nonanalytic for |c| — 1. 
This is also visible from momentum space perspective 
through discontinuity of /fe at /cc = and |c| = 1. Indeed, 
it is easy to check that 

linife_^fc, /fe = -1 for |c| < 1, 
limfc_>fe^ /fc = for |c| = 1, 
limfe^fe^ /fc = 1 for |c| > 1. 

This discontinuity is presented in detail in Fig. [TOl 

Appearance of oscillations and some constant prefac- 
tors in ((65)) . overlooked by "continuous" J-"smooth, can be 
explained in a similar way as in Sees. IIII Al and IIIIBI 
Again they can be traced back to the fact that crosses 
zero when |c| < 1 (Fig. [T5|) . This leads to logarithmic 
divergence of In causing corrections to J^3mooth• 

More quantitatively, when |c| = 1 then reaches zero 
for /c = /cc = 0. Similarly as in Sec. IIII Al it leads to 
the additional term, saturating at In2/2A^, which has to 
be added to the right-hand side of ([75)1 . It results in the 
prefactor of \/2 in (p5|) . When |c| < 1, /fe reaches zero 
for fco ^ 0, which is different than fee = for which the 
gap closes. It is easy to check that 

^« = ™ (irlJIrr^) " 2i<5i/r^. (74) 

This fco is in general incommensurate with momenta 
k ([68]) and it leads to oscillations of fidelity. Wc 
can directly use the formalism presented in Sec. IIII Bl 
to obtain the additional factor of 2 |cos(fco-^/2)| ~ 
2 |cos(|(5|A^\/l — c^)! modifying J'smooth- Note that this 
is the prefactor appearing in (|65p and that the approxi- 
mation (|74p has been used to simplify cos(fcoA^/2). 



V. FIDELITY AND NONEQUILIBRIUM 
QUANTUM QUENCHES 

Our goal here is to briefly illustrate how our results 
and analytical techniques presented in the previous sec- 
tions can contribute to a better theoretical description 
of quantum quenches (both instantaneous and continues 
ones). 

The simplest quench one can consider is the instanta- 
neous quench, where the parameter of the Hamiltonian, 
say A, is changed at once from Ai to A2. Assuming that 
the system was initially prepared in the ground state | Ai) , 
the modulus of the overlap of its wave-function onto the 
new ground state IA2) reads (P) 

|(Ai|A2)| = J-(Ai/2 + A2/2, A2/2 - Ai/2). (75) 

In other words, fidelity provides here a square root of 
the probability of finding the system in the new ground 
state after the instantaneous quench. If the system is in 
the thermodynamic limit and (5 = Ai/2 — A2/2 is small 
enough, depending on A = Ai/2 -f A2/2 either ^ or 
(O can be used to evaluate ([75t near the critical point 
(A« A,). 

Another application of the analytical techniques de- 
veloped in the former sections comes in the context of 
the density of quasipa rticles excited during an instan- 
taneous quench pH - [l6{ . For clarity of our discussion we 
will focus here on the instantaneous transition in the Ising 
model, where the system is suddenly moved from (51,71) 
to (52,72) with (7i_2 = 1 -I- e ± (5 and 71,2 = 1 (path A 
on Fig. [TJ see ^ for the Hamiltonian). The density of 
quasiparticles reads 

- r - vT = r -(1 - fi) = mic) + o{n 

Jo Jo TT 

(76) 

where p'f^ is the probability of excitation of the A:-th mo- 
mentum mode during the quench [TgI . l43j , fk is related 
to fidelity ([75]) through ([ MT^ . and again we assume \S\, 
|e| < 1 and c = e/\S\. 

The scaling function B{c) - presented in Fig. [^Dk - is 
given by 

^ (1 - \c\)lmE{c2) + 2K{c,) for |c| < 1, 

2tt \ (|c| - l)Imi;(c2) - 2K{ci) for |c| > 1, 

(77) 

where ci, C2, K and E are given by ([^ and (15(1)) . To 
obtain this result we used the approximation p5|). 
Equation d7S|) extends the results of [ll| where density 
of excited quasiparticles Uex was calculated for a specific 
situation where the instantaneous quench has either be- 
gan from the critical point or ended at the critical point. 
Note that we consider the instantaneous quench every- 
where around the critical point. Moreover, our calcula- 
tion provides the scaling function B(c), which is indis- 
pensable for capturing full universal description of n^x 
for any two states close to the critical point. It also al- 
lows for unified description of the density of quasiparti- 
cles for instantaneous quenches both near and away from 
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FIG. 20: (Color online). Scaling function (|77|) . and its deriva- 
tive (|78p . characterizing density of quasiparticles after an in- 
stantaneous quench in the Ising model (I76p . Panel (a): the 
black line shows analytical expression (|77|) . while red crosses 
show numerically evaluated Wex/I^l from (|76|) for TV — 10^ 
and 5 = 10~^. Panel (b): logarithmic divergence oi dB{c)/dc 
around c — 1 (|78p . The solid black line shows analytical result 
(|78|l . The dashed lines present numerics done for 5 = 10""^ 
and N equal to 10^ (red), 2 x 10^ (blue) and 4 x 10^ (green) 
[top to bottom] . The logarithmic singularity is reached in the 
limit of A'' — >■ oo (compare with Fig. [4] and discussion in Sec. 

MM- 



the critical point. To illustrate the latter, we consider 
|c| > 1, but still < |e| < 1, and expand ^ getting 
B{c) ~ l/4|c|. Putting it into ([751) we obtain 

Hex ^ S^/A\e\. 

We also note that the scaling function B{c) is nonanalytic 
for |c| — 1. This results from discontinuity of fk, and 
consequently p™, at fc^ = and |c| = 1 ([55)) . Indeed, the 
derivative of B(c) can be expanded around c — 1 



dB{c) 



dc 



2-31n2 ln|l-c| 



2tt 



27r 



+ 0{c-l). (78) 



This logarithmic divergence is reached when the size of 
the system iV — >■ oo (see Fig. [^Ob and the related discus- 
sion in Sec. IIIIAp . 

A more complicated quench dynamics shows when the 
parameter driving the transition, say A, is gradually 
shifted 

X{t) = Ac + t/TQ, t = —oo . . . CO, 

where the quench time scale tq controls how fast the 
system is driven: see [l^ for a review of the resulting 
quantum dynamics from the perspective relevant for our 
discussion. Briefly, it turns out that the system's evo- 
lution can be approximately divided into three regimes 



adiabatic 
direction of evolution 



impulse 



adiabatic 



X -X 



x+x 



FIG. 21: Schematic illustration of dynamics of a generic quan- 
tum phase transition (see text for details). 



(Fig. [?!]) fij, |4^- First, the adiabatic regime takes place 
as long as the system is far enough from the critical point: 
A(t) < Ac — A, where A is a small parameter that will be 
discussed below. Assuming that the evolution started 
from a ground state, wave-function of the driven system 
is given here by the instantaneous ground state of its 
Hamiltonian. Then the impulse (diabatic) regime hap- 
pens around the critical point, X{t) G (Ac — A, Ac + A), and 
system's wave-function is assumed not to change here. 
Therefore, it is the same as at the last "adiabatic" in- 
stant t such that X{—t) — Ac — A. Finally, the system 
enters the second adiabatic regime when A(t) > Ac + A. 
It is assumed that no additional excitations are created 
by the quench in this regime. The parameter A is pro- 
vided by the quantum version of the Kibble-Zurck theory 

M 



A 



where tq ^ 1 is assumed (the slow quench limit) and 
the prefactor omitted above is 0(1): A ^ 1. In addi- 
tion, one also assumes here that the correlation length 
(^) during the quench is much smaller than the system 
size to avoid finite size effects. In the adiabatic-impulse 
approximation this implies that L ^ ^(Ac — A) ~ X~'^ or 
simply N/tq^^^^^'"^ > 1, where N = L'^. 

It turns out that one can relate all these concepts to 
fidelity and use ^ to quantify scaling of the system exci- 
tation with the quench rate tq . Employing the adiabatic- 
impulse approximation, we see that the population left 
in the instantaneous ground state away from the critical 
point - i.e., after leaving the impulse regime - is given 
by 



|(Ac-A|Ac + A)|2=^2(^^^^)^ g 
or, after using (HI), by 

(^~N X const/rg 



-l/(l+zi^) 



exp I 



d!//(l+ziy) 



(79) 



where the const prefactor is given by 0(2^(0)). This 
result has a simple interpretation in the context of sym- 
metry breaking phase transitions. As predicted by the 
Kibble-Zurek theory, such transitions lead to creation of 
topological defects whose density is given by 
[2^ . Thus, the probability of finding the system after 
the quench in its instantaneous ground state is exponen- 
tial in the number of topological defects created during 
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the quench (see Sec. V of [46| for the same conclusion 
worked out in the Ising model; note that ([7^ is much 
smaller than unity when the adiabatic-impulse approxi- 
mation applies). 

This is in agreement with studies of dynamics of the 
quantum Ising chain, whose Hamiltonian is given by ^ 
with 7 = 1 (note that g in ([6]) corresponds to A here) 
[38.. J:6]. In Sec. V of |46i] the Ising chain was initially 
exposed to a large magnetic field: evolution begun from 
a paramagnetic phase ground state. The magnetic field 
was ramped down on a time scale tq and the system 
was moved across the paramagnetic-ferromagnetic criti- 
cal point. The evolution stopped at the zero magnetic 
field (deeply in the second adiabatic regime) and the 
squared overlap between the evolved wave-function and 
the ground state at the zero magnetic field was found to 
equal 

exp (— iV X const/y^) , 

where const = 0(1). This is in agreement with ((79)) 
because = z = d = 1 for the quenched Ising model. 
Naturally, it would be also interesting to verify our simple 
scaling result (l79t in other models as well. 
At the risk of laboring the obvious, we note that fidelity - 
or more interestingly the scaling factor d pi[) - is insen- 
sitive to both the relaxation processes after the instanta- 
neous quench and the adiabatic dynamics following cross- 
ing of the critical point. This shall make it a useful tool 
for theoretical characterization of quantum quenches. 



VI. CONCLUSIONS 

We have extensively studied ground state fidelity in the 
XY model illustrating its rich behavior around different 
critical lines. This supports its potential applications as 
an insightful probe of quantum criticality 8j. A special 
stress has been placed on discussion of fidelity from the 
momentum space angle allowing for better understanding 
of the studied models. We have derived a variety of new 
analytical results in the thermodynamic limit and verified 
them through numerical simulations. 

First, we have discussed in detail an example where 
fidelity is dominated by the universal contribution cap- 
tured by the scaling relations (|4l5p derived in [l3| (Sec. 
nil Ap . Second, we have focused on the case, where close 
to the critical point the "smooth" universal contribution 
to fidelity is strongly modulated: it is multiplied by an 
oscillating factor (Sec. IIIIBI) . This interesting effect has 
been qualitatively explained and accurately analytically 
described. Third, we have characterized fidelity calcu- 
lated along one of the critical lines where non-universal 
contributions provide the key input (Sec. IIIICp . Fourth, 
we have discussed fidelity near the multicritical point, 
where the universal scaling relations are violated and a 
more advanced scaling theory needs to be worked out 
(Sec. ImPl) . 



The techniques used in this article may be applicable 
to other systems as well. They have been tested out- 
side of the XY model in the extended Ising model (Sec. 
lIVp . They have also been applied to get some additional 
analytical insight into dynamics of quantum phase tran- 
sitions (Sec. N} . 

From the technical aspects of our work, we would like 
to stress non-triviality of the exchange of the product 
of over momentum contributions into an integral over 
their logarithms. More fundamentally, we mention non- 
commutativity oi N ^ oo and (5 — >■ limits near the 
critical point. Depending on the order in which they are 
taken, one ends in either the "small system" limit or the 
thermodynamic limit and fidelity behaves very differently 
in these two limits around the critical points. 
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Appendix : Estimation of the errors of analytical 
approximations 

Calculation of the leading universal contribution to fi- 
delity in Sees. IIII Al and llll BI reauires integration of In fk 
over k from to tt. Even though /^'s are known exactly 
(llimsp . such an integral cannot be analytically done to 
the best of our knowledge. To proceed, we have approx- 
imated In/fc in Sees. UTTH and HITB] bv (^5]) and P^ . 
respectively, and extended the integration range over k. 
Below we estimate the errors resulting from these ap- 
proximations. 

Approximations in Sec. IIII Ai We demonstrate 
here that the difference between our analytical expression 
(l27l) and the exact result scales as 0{6'^/j^). Namely, we 
are going to show that 



2tt 



dkiiifk~{-mic)/^) 



(Al) 



is 0(^^/7^). The first term on the right-hand side of (lAip 
is the exact expression [fk is given there by (|lll23l24p ]. 
while the second term is its approximation (P7|) . Note 
that we do not study here the errors coming from ex- 
change of the product over fk into the integral (ITSt . 

For simplicity, we restrict our analytical studies to the 
specific case when one of the states is exactly at the crit- 
ical point: e = (5 > (c = 1). Straightforward, but 
tedious, generalization of the proof listed below allows to 
extend it to any |e| ^ 7^. A numerical check that the 
same error appears when e 7^ (5 (c 7^ 1) is illustrated in 

Fig. m 

Error (jAip is bounded: \E\ < Ei + E2, where 



El 



2-K 



dk In fk 



2-K 



dk In fk 



< — dk 
- 27r 



In fk - In fk 



(A2) 
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FIG. 22: (Color online). Numerical calculation of the error 
(|A1[) for different 7, 5 and c = t/\5\. All the curves show 
that \E\ X 7^/(5^ is bounded from above by about 0.2. They 
also illustrate that E scales as /■y'^. Upper panel: S = 10"'^ 
(dashed red line) and S — 10~^ (solid black line); 7 = 1 for 
both curves. Middle panel: 5 = 10"* (dashed red line) and 
5 = 10~^ (solid black line); 7 = 0.5 for both curves. Bottom 
panel: 5 = 10^^ (dashed red line) and 5 = 10~^ (solid black 
line); 7 — 0.1 for both curves. We keep |e|, \S\ ^ 7^ in these 
calculations. 



is due to approximation of pk and qk in ((25)) by pk and 
Qk (see below) and 





1 f°° 


1 




E2 = 


— / dklnfk 


< — dk 


In A 




27r 


~ 27r^ 



results from sending the upper integration limit in (|26|) 
to infinity. For clarity of our discussion, we repeat the 
expressions for exact fk 

Pk = (1 - cosA;)(l - cos/c + 2(5) + 7^ sin^ A:, 
qk = 2(57 sin fc, 

and approximated fk 



inA ^ im i 



Pk 



2v^l+l, 



In order to estimate Ei we find a function such that 



In fk - In 



(A3) 



To further simplify our calculations, we assume that 
(5/7^ < 0.2: other bounds would give us a slightly dif- 
ferent pref actor in the estimation of the magnitude of E, 



but will not change its scaling with 7 and 5. To derive 
(|A3p we start with inequality: 



ln/fe-ln/fe| <0.6 AV/I-1 



(A4) 



which is true when fk/ fk ~ 
proved below). We consider 

fk/ Ik ~ 



< 0.25 (the latter wih be 



1 



(A5) 



Two remarks are in order now. First, we will use below 
that Pk,qk,Pk,qk > and < Pk/qk < 00. Second, we 
will consider two regions in k separately: 

• k e [0, 2(5/7] or equivalently < Pk/qk < 1- 
In this case we estimate (IA5I) as 



fk/.fi~i 



Pk/qk 



Pk/qk 



1 



Pk/qk 



< 



Pk/qk Pk/qk 



< \Pk/qk - Pk/qk\ = .gfc/0.6. 



Above, in the first step we bounded denominator 
by 1 and in the second step we used inequality 

x/VT 



- y/ 



'I 



y 



< \x ~ y\. Subsequent study 



< 0.25 when 



of \pk/qk - Pk/qk\ reveals that f^/f^ - 1 

(5/7^ < 0.2. This in turn validates inequality (|A4p in 
the considered fc-region. The above estimation provides 
us also with the expression for gk from (IA3|) : gk — 
0.6\pk/qk —Pk/qk\- We integrate it over k, and simplify 
the result getting 



25/7 

/ dkgk <o.n ^Vt' 



1 

2^ 



• fee [2(5/7, or equivalently < qk/pk < 1- 
In that case we estimate (IA5I) as 



(A6) 



fk/fk 1 



1 



< 



< 



1 



1 + 1/V2 
1 



1 



V2 



ql/Pk - fk/pl\ 



5fc/0.6, 



where in the first step we bounded denomina- 
tor by 1 + I/V2 and later we used inequality 

2 



i/VT 



r 



< 



y \/2. A careful 



19 



0.25 



to 



0.248 



0.246 




Estimation of E2 is straightforward because 



E,,< — 




1 



2 ■ 

52 



47r7 J^^/s 47r27' 



< 0.03 (577^ 



(A9) 



< 



Summing up the contributions (|A8|) and ()A9p we find 
that 



FIG. 23: Integration error (|A10|) . The curves from bottom to 
top correspond to 5 = 10~^, 10~* and 10~^, respectively. For 
each 5 we choose g = 0,0.9, 1. The curves for different g's, 
but the same 5, perfectly collapse onto each other. This plot 
supports E ^ 0.25(5^. 



study of 



2+V2 



\iI/pI - €IpI\ shows that fl/fl - 1 



< 



0.25 when 5/7^ < 0.2. It vahdates inequality (|A4I) in 
the studied /c-region. From the above calculation we see 



that gk from (|A3p is equal to 



0.6 



Tk/Pk - TklP. 



fe- 



lt 

can be integrated over k and the obtained result can be 
simplified to yield 

^ r dk gk < 0.25 Sy-f^. (A7) 



/ 2.5/7 

This concludes our proof of 
Combining with dM 



El < 0.36 (577^- 



a for 5/7^ < 0.2. 

H and (IA7I) we obtain 



(A8) 



\E\ < 0.4(577^ 



when e = (5 (c = 1). This is in agreement with Fig. [521 
which in fact suggests a stronger bound, and extends our 
analytical results beyond the e — S case. 

Approximations in Sec. IIII Bl We justify here 
numerically the 0{S^) error in To that end we look 
at 



E=^ J\klnfk- i-mMc)) , (AlO) 

where this time fk is given exactly by (|lll40l4ip . The 
numerical results for various values of parameters are 
presented in Fig. [23l It shows that the error scales as 
0{6^) confirming that our approximation correctly cap- 
tures leading universal contribution to J^3mooth■ We note 
also that our numerics suggests that the error, and in 
principle the exact value of the integral (|44|) . do not de- 
pend on 5 S [—1,1]. 
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